Boundary condition for Ginzburg-Landau theory of superconducting layers 
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Electrostatic charging changes the critical temperature of superconducting thin layers. To un- 
derstand the basic mechanism, it is possible to use the Ginzburg-Landau theory with the boundary 
condition derived by de Gennes from the BCS theory. Here we show that a similar boundary con- 
dition can be obtained from the principle of minimum free energy. We compare the two boundary 
conditions and use the Budd-Vannimenus theorem as a test of approximations. 

PACS numbers: 74.20.De 



INTRODUCTION 



Much experimental effort is devoted to find supercon- 
ducting materials with critical temperatures as high as 
possible. It is well known that the critical temperature 
depends on the charge carrier density. The charge car- 
rier density can be changed by doping and to some extend 
it can also be changed by electrostatic charging. Conse- 
quently it is an attractive task to determine, how electro- 
static charging evoked by an applied electric field changes 
the critical temperature of superconductors. Generally 
the experiments revealed that it is easier to increase T c 
than to decrease The Ginzburg-Landau (GL) the- 

ory with the de Gennes boundary condition can be used 
to understand this behavior—. 

The aim of this paper is to show how a superconductor 
screens the external electric field and to which extent the 
boundary condition derived from the minimum free en- 
ergy principle is compatible with the de Gennes boundary 
condition. 



Charges at a solid surface partially leak out of the 
surface. This creates a surface dipole. The Budd- 
Vannimenus theorem^ describes the step in the surface 
potential due to this surface dipole as a simple expres- 
sion of the bulk free energy density. Therefore it is well 
suited to test the approximations used in this paper. 

In the first chapter we explain the model and the parts 
considered in the free energy of the superconductor and 
solve the Euler-Lagrange equations for the GL and charge 
carrier wave function and the surface potential. In chap- 
ter III the corresponding equations outside the super- 
conductor are solved and in chapter IV the continuity 
requirements determine the remaining constants. Chap- 
ter V presents the numerical values which are compared 
with the de Gennes boundary condition in chapter VI. 
Finally we conclude in chapter VII. 



II. FREE ENERGY IN THE 
SUPERCONDUCTOR 

We start with the free energy 

T = J (/tf + fee + /gl + /cist) dr, (1) 

where we include only the terms most relevant for the 
above specified problem. 

The first term /tf is the Thomas-Fermi internal en- 
ergy, for which we use the LDA (local-density approxi- 
mation) 



/tf = r (3tt 2 ) 3 — ri3. 
5 2m 



(2) 



The second term /qc represents the condensation energy 
for which we use the formula following from the Gorter- 
Casimir two fluid model 5 -^ 
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The electrostatic energy density term reads 



/cist — 2 e °^ 2 



eipSn, 



(3) 



(4) 



in the form suitable for performing variations. For sim- 
plicity we exclude the magnetic field and it's related ki- 
netic energy of the screening current. We take the vector 
potential A to be zero and write the GL gradient term 
as 



/ GL = A^| W | 2 



(•5) 



Here we have chosen the GL wave function if) normalized 
with respect to the total charge carrier density n. In 
the spirit of the Thomas- Fermi approximation the charge 
carriers are described by a wave function vb n with n = xjj 2 
and the superconducting fluid density used in the formula 
for the condensation energy (|3|) reads 
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In short, the free energy is expressed by three inde- 
pendent variables: the scalar potential p determining the 
electric field E = — Vy>, the GL wave function ip and the 
charge carrier wave function tp n . We assume that the 
material parameters, the critical temperature T c and the 
Sommcrfeld parameter 7 depend on the charge carrier's 
density n, by using the approximations 

T c (n) = T c0 fl + ^i§^ 
Tliat ohm 



(8) 



where m a t is the crystal lattice density. In the following 
we shall use three characteristic length: 

(i) the Thomas-Fermi screening length A^ F = 7^%Ep, 

(ii) the Bohr radius Or = ine °!} 

(iii) and the coherence length En = 1 — ^a at ■ 

From the charge neutrality requirement we know that 
^noo = y/n\at (here and in the following the subscript 00 
denotes the magnitude far from the surface). To keep 
things simple, we use the approximations 



"0n = (1 + Slp n ) 

V> = Voo(l+<V) 
eip = E F (ifao + Sip) 



(9) 
(10) 
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and suppose that the deviations of the three indepen- 
dent variables from the optimum values are small. In 
a homogeneous superconductor far from the surface all 
these deviations have a zero value and the derivatives of 
the free energy formula ([1]) with respect to them must 
be also zero. From these requirements we get the magni- 
tudes of the optimum superfluid fraction 
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SiJj = 0, (16) 
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Sip = 0. (17) 
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Close to the planar surface we can assume exponential 
dependencies of the deviations and from the EL equations 
(|15H17p we get a second order equation for the square of 
the expected penetration depth A. Two solutions arise 
out of it. 

Observing that Atf^b <C £,q, we find a first approxi- 
mate solution A in the form of the coherence-like length 



6 = £o 
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(18) 



In this solution the scalar potential is constant and local 
charge neutrality is preserved. Only the deviation of the 
wave function Sip is nonzero (Cj will denote its magni- 
tude). 

A very small penetration depth characterizes the sec- 
ond solution such that this solution can be simplified. 
Using the same approximation as above, we find that 
the second penetration depth equals the Thomas-Fermi 
screening length Atf- In this solution the scalar poten- 
tial displays a sharp step (Ctf will denote its magni- 
tude) and from the Poisson equation (TTSJ) follows that the 
charge carrier's density changes accordingly. The sharp 
step on the GL wave function ip is negligibly small due 
to the factor A TF /£ 2 which enters the resulting formula. 
It corresponds to the well known fact that the GL wave 
function ip cannot abruptly change. 

The general solution can be written as a sum of the 
two above described solutions: 



and the optimum magnitude of the scalar potential 
(_ „s dhxT 



Poo = -1- 



r2 3§2(2(l-i 4 )^ + (l + t 



dlnn 



<91nn / 
(13) 

The electrostatic potential energy of the charge carrier 
thus equals the Fermi energy 



ft 2 , 2 ,2/3 
Ei-p = — o7r n) 
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(14) 



with a small (lower than gap) correction represented by 
the second term in (fTB"]) . 

Using the second order expansion of the free energy 
([1]) , from the variation we get three linear Euler-Lagrange 
(EL) equations for the three independent variables, 
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Here Ctf describes the step of the scalar potential in 
units of Ep/e according to (fTTj) . Using (TT9")) - (|2"Tj) we can 
calculate the free energy 
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For a semi-infinite medium the first term gives an infinite 
contribution which is not influenced by the surface con- 
ditions, so we do not need to deal with it. The last three 
terms correspond to the surface energy, which according 
to the principle of minimum free energy should take an 
extremum. The minimum of the free energy is obtained 
for 



3A T f6 91nT c 



dhin 
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(23) 



in which case the derivative of the wave function ip at 
the surface is zero. For lead at t = 0.9 we get = 
— 0.00044Ctf- As expected, the deviation of the GL wave 
function Sip is much smaller compared to the sharp steps 
on the scalar potential Sip and on the charge carrier wave 
function Sip n . 

We see that the principle of minimum free energy en- 
tails the GL boundary condition. Towards the surface 
the GL wave function Sip displays a small gradual change, 
only very close to the surface its derivative jumps to zero. 
The solution is complete, if the parameter Ctf is deter- 
mined, ft can be derived from the requirement of con- 
tinuity with a solution minimizing the total free energy 
including the one of the vacuum outside. 



III. FREE ENERGY OUTSIDE THE SLAB 

Now we approximate the free energy density outside 
the superconductor by 



F = / (/w + Jgl + /cist) dr. 



(24) 



We include the electrostatic term, the GL gradient cor- 
rection and the von Weizsacker kinetic energy functional 



/w 



2 m 



(25) 



fn the limit of rapidly varying densities this kinetic en- 
ergy term is dominant and when describing charge carri- 
ers tunneling outside the superconductor this term can- 
not be omitted. We have not included this term into the 
formula ([5]) describing the free energy inside. The reason 
is that inside the superconductor the Thomas-Fermi in- 
ternal energy plays the dominant role and moreover, as 
it is shown e.g. in the book of Dreizler and Gross^, in 
the limit of nearly homogeneous systems the second order 
term of the gradient expansion provides a better approx- 
imation. It has the same structure as the von Weizsacker 
kinetic energy functional, but its coefficient is nine times 
lower. We suppose that for the rough estimates presented 
here this relatively small correction can be neglected. 

In the vacuum far from the surface the scalar potential 
reaches the magnitude of the work function <pw, so that 
we can approximate 



eip = E F ((fiw + S<p) 



(26) 



The density of the tunneling charge carriers quickly drops 
to zero. Using an analogous notation as above we write 



(27) 



supposedly that Sip n is small. For the wave function ip 
we use the approximation 



1p = Ipoc + Slp, 



(28) 



where ipoo represents the supcrfiuid fraction in the vac- 
uum far from the surface. Let us remind that ip is nor- 
malized to the charge carrier density, see ((6]), so that 
does not need to be zero. The free energy density in the 
vacuum thus reads 



/out = ^Pf (2 (V<5V„) 2 + Hi W 



-A^ F (VScp) 2 + {<p w + S<p) Sijjl 



(29) 



and we can write the Euler-Lagrange equations. 

The variation with respect to the wave function ip gives 
the condition 



S^lV 2 8iP = 0. 



(30) 



We see that Sip remains constant or changes linearly. 

The proximity effects indicate that the correlated 
charge carriers can remain correlated even if they are 
tunneling. For simplicity we suppose, that the superfluid 
fraction of the charge carriers tunneling outside the ma- 
terial does not change, so we take Sip = ip(Q). The two 
other Euler-Lagrange equations read 



2(fi W Slp n 
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W 2 Si/j n = 



Hi + g A TF V 2 ^ = 0. 



(31) 



(32) 



In the same way as above we can try the exponential 
solution 



6ip n = K n exp 
Sip = K v exp 



- 

Aw/ 

2x 
Aw 



(33) 
(34) 



where Aw denotes the tunneling length which follows 
from the Euler-Lagrange equation (f3"2"| as 



A 



w 



K n 



The work function can be determined from (1311) as 



I6A^ F 
7r 2 a!A w ' 



(35) 



(36) 



In this way we have an approximate solution outside the 
superconductor, which should be linked to the solution 
inside. 



IV. CONTINUITY REQUIREMENTS 

At the surface the continuity of the wave function tp n 
and the continuity of the scalar potential tp with its 
derivative (continuity of the electric field) must be en- 
sured. We get three conditions 
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(37) 



(38) 



(39) 



where the term E a representing the applied electric field 
is included into the condition of continuity for the electric 
field. From the continuity requirements (|37H39[) we obtain 
the equation 
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1447r 2 a| (Ctf - E a 



(-4 + 3C TF ) 
-1 + Ctf = 0, 



(40) 



determining the step of the scalar potential Ctf ■ 



NUMERICAL VALUES 



The sixth order equation (I40|) can be numerically 
solved. For small applied electric fields the linear ex- 
pansion 



Ctf — C' 



TFO 



(41) 



is applicable and for lead at temperature t=0.9 we get 
the numerical solution 



C TF = 0.457 - 0.53£ a 



(42) 



The numerical estimate for the tunneling length follows 
to be from (|35|) Aw = 3.17 Atf and the work function 
according to (|36|) pw— 1-43 eV. Taking into account how 
many simplifications we have used, it is surprising that 
the obtained results seem to be quite reasonable. The 
estimated magnitude of the work function is comparable 
with the experimentally determined value of pw= 4.25 
eV4 

The sharp step of the scalar potential can be estimated 
from the modified Budd-Vannimenus theorem^ according 
to which 



Here f c \ denotes the spatial density of the electronic 
free energy, which can be roughly approximated by the 
Thomas- Fermi internal energy /tf defined in (J2j) . Then 
the Budd-Vannimenus theorem (|43[) predicts a sharp 
step, Ctf = §> of the scalar potential in units of Ep/e. 



The numerical solution (|42|) gives a comparable result 
what strongly supports the applicability of the here used 
approximations. 

We saw that the numerical values of the measurable 
quantities are reasonable. In Fig. [1] the scalar po- 
tential is plotted. As expected, inside the supercon- 
ductor the scalar potential acquires the Fermi energy 
value, while in the vacuum outside it reaches the work 
function value tp w . The dashed and dotted lines corre- 
spond to the experimentally accessible applied electric 
field E a = ±0.01^ S ±1.7 x 10 7 V/cm. As it is seen 
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FIG. 1: The electrostatic potential tp calculated for lead at 
t = 0.9K. The dashed and dotted lines correspond to applied 
electric field E a = ±0.01 

in Fig. [21 the external electric field is screened on the 
Thomas-Fermi screening length. In the figure [3] devia- 
tions of the charge carrier densities from the equilibrium 
values are plotted. We can see that the magnitudes of 
these deviations are small. Close to the surface the su- 
pcrfluid density n s decreases and this decrease is com- 
pensated by an increase of the normal fluid density n n . 
The total charge carrier density n shows no change on 
the scale of the coherence length. 



VI. COMPARISONS WITH DE GENNES 
FORMULA 

Now we compare the GL boundary condition follow- 
ing from the minimum free energy principle with the de 
Gennes boundary condition^ 
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FIG. 2: Charge density p. To make the screening visible, the 
dashed and dotted lines correspond to the applied electric 
field E a = ±0.3 
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FIG. 3: Deviations of charge carrier densities. To ensure vis- 
ibility, the dashed and dotted lines correspond to the applied 
field E a = ±0.1. 



according to which the derivative of the gap at the surface 
is not exactly zero even without external electric field. 
The zero field extrapolation length bo is around 1 cm 
(almost infinity from the microscopic point of view). The 
effective potential U s 



1 
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(45) 



determines how the extrapolation length b changes if an 
external electric field E a is applied 3 . De Gennes esti- 
mated the surface ratio 



A(0) 
Ac 



(46) 
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FIG. 4: Deviation Sip of the wave function if) from the equi- 
librium value. Close to the surface its derivative approaches 
zero (see the insert). The dashed and dotted lines indicate 
how the extrapolation length changes with applied field. 



to be close to one. For lead the formula (|45|) gives U s = 
1.35 x 10 7 V. From the minimum free energy we know, 
however, that the derivative at the surface should be zero. 
In Fig[4]we see how the deviation of the wave function ip 
at the surface decreases with the derivative determined 
by the parameter Cj. From this we get the extrapolation 
length bo pa 2.8 mm, a value comparable with the one 
estimated by de Gennes. Only very close to the surface 
(on the distance of Thomas-Fermi screening length) the 
derivative of the wave function ip approaches zero (see 
insert of the Fig. [4j. 

In figured] we can also observe how the extrapolation 
length changes if an electric field E a is applied. By sub- 
stituting b — — ^t/Cj into (|4"4"|) and using (f2"3"| with the 
approximation (|41[) we get a simple expression for the 
effective potential U s 
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This formula is similar to the de Gennes formula (|45|) . 
We should notice, however, that in this formula the ex- 
trapolation parameter £/2 of (T4"Tj) appears instead of the 
surface ratio i] which enters de Gennes formula (|45|) . 



VII. CONCLUSIONS 

It was shown in this paper that the minimum free en- 
ergy principle entails a zero derivative of the wave func- 
tion ip at the surface of the superconductor. On the 
scale of the coherence length, however, even if no ex- 
ternal electric field is applied, the derivative is nonzero 



6 



and its magnitude corresponds to the de Gennes esti- 
mate. Only on the Thomas-Fermi screening length scale 
it approaches zero. In the presence of an external elec- 
tric field the extrapolation length changes according to 
equation (J44J) , with the effective potential given by equa- 
tion (|47|) . This formula is similar to formula (|45| follow- 
ing from the de Gennes theory. The agreement with the 
Budd-Vannimenus theorem and the numerical estimates 
support the applicability of the proposed approach. 
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